A NEW FORMALISM FOR THE COMPUTATION OF RKKY INTERACTION IN APERIODIC SYSTEMS 
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A numerical method to investigate RKKY interaction between localized spins in aperiodic materials 
is presented. Based on an expansion of the spectral measure on orthogonal polynomials, this ap- 
proach leads to an effective framework to analyse linear response formula for the RKKY coupling in 
complex systems. Particularly useful in the tight-binding scheme it is used in this paper to probe the 
role of the local environment and the quasiperiodic potential on the interaction between magnetic 
spins. Interesting features are revealed and discussed within the context of anomalous localization 
and transport. 



PACS numbers: 72.90.+y 61.44.Br 72.10.-d 
I. INTRODUCTION 



The purpose of this work is to present a new for- 
malism based on real space recursion scheme |l| which 
enables to compute the so-called Rudermann Kittel Ka- 
suya Yosida interaction (RKKY) [Q-g an effective cou- 
pling between two localized magnetic moments, mediated 
by conduction electrons. 

This long range oscillatory interaction is now well 
understood in pure metals and in weakly disordered sys- 
tems, or even close to the Anderson transition. RKKY in- 
teraction has been the subjected of a great attention dur- 
ing the past years for instance for understanding the spin 
glass transition H, magnetic long range order in high- 
T c cuprates M, or more recently because of its anoma- 
lous behavior related to giant magnetoresistance effects 
in magnetic multilayers S. 

As conduction electrons carry this interaction, one 
can wonder what happens for systems in which electronic 
propagation is anomalous B or in cases of inhomoge- 
neous disorder ? (e.g. local inhomogeneities, correlated 
disorder,...). These issues may be addressed in relation 
to unusual RKKY interaction in quasicrystals, where the 
amplitude of the coupling was found to be anomalously 
strong and with similar values for several manganese- 
based quasicrystals [jllj with different densities of states. 
Such results were found to apparently contradict the 
classical theory stating that RKKY should be propor- 
tional to the total DOS at Fermi level. 

In the following, after reminding some general fea- 
tures of RKKY coupling, the method will be detailed 
followed by a first part addressing the numerical conver- 
gence of the algorithm. A second part is devoted to the 
study of the impact of local disorder or quasiperiodic po- 
tential on RKKY. This will exemplify the interest of the 
method. 



II. CALCULATION OF THE NON-LOCAL 
SUSCEPTIBILITY 



Indirect RKKY interaction stems from the coupling 
between localized magnetic moments and propagating 
electrons. If an electron in a state of energy E < Ep 
undergoes a transition to a state of energy E' > Ep be- 
cause of the coupling with the localized moment in | n), 
then a hole-electron pair is created and propagate coher- 
ently during a certain time r, with \E' — E\ < Ti/t, until 
the pair is destroyed by diffusion on another magnetic 
impurity located in | rj). Consequently the longer is the 
propagation time the smallest will be the vicinity around 
Fermi energy that account for RKKY. The generic form 
of the effective coupling between two magnetic impurities 
mediated by itinerant electrons reads : 



Zrkky (n , rj , E) = J 2 x(ri , r,- , E)S Ti -S rj 

with J is the interaction between the localized mo- 
ment S ri and the spin of the itinerant electrons, and 
x(ri,rj, Ep) contains the sum of all the electron-hole 
propagation paths from \ri) to \rj). RKKY is then pro- 
portional to the electronic susceptibility x( r ij) of itiner- 
ant electrons. When J > (resp. J < 0), the configura- 
tions of parallel spins (resp. antiparallel) will minimize 
the energy promoting ferromagnetic state (resp. antifcr- 
romagnetic state). The susceptibility as a contribution 
of all the scattering pathes of the hole-electron pair can 
be written down as 



x{n,rj 
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introducing the retarded (G_) and advanced Green's 
functions (G+) which define the amplitude of propaga- 
tion of the hole-electron pair. Note that there exists 
an exact sum rule between the susceptibility and the 
local density of states (LDoS) J2jX( r i, r j) = Pi{E) = 



^m( n I G(z) 
test. 



which could be used as a numerical 



In metallic systems with space dimension D, the 
interaction is given by 

1rkky{t, E f ) ~ A{t) cos(2k F r + 5{v))/r D 



(1) 



which manifests a long range oscillating behavior. pf For 
a free electron gas, the A, 8 are independent of r, whereas 
for weak disorder limit, A(r) becomes a random but 
smooth function of r and <5(r) is the phase shift associated 
to the scattering of electrons on impurities. <5(r) becomes 
random for r > l pm (mean free path) . It has been shown 
however that the mean free path is not a good measure 
of the typical length scale of this interaction. Indeed, 



the disorder average of even moments (x p (l 



I)) 



will contribute roughly as (\/r D ) 1 ' 2p , whereas odd mo- 
ments of the susceptibility will vanish exponentially over 
the mean free path length as expressed by the following 
averaged results [|| 
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{X p+1 (\r t -r,\)) ~ exp(- 



taking p(Ep) the DoS at Fermi level, l m the mean free 
path, and fi p a constant independent of the parameter of 
the hamiltonian. Consequently, the average over disorder 
of electronic susceptibility is not sufficient for describing 
the correct range, amplitude and phase of the interaction, 
which fluctuates very much from one random configura- 
tion to another. MM 

Let us now consider a simple argument for the 
RKKY coupling. If one considers a cube of length L 
(with periodic boundary conditions) which is larger than 
the typical distance magnetic impurities, then the expan- 
sion of the susceptibility in a basis of eigenstates {|n)}, 
reads : 



X (r i ,r j ) = 2^eJ2 



(ri\n)(n\ Tj )(r j \m)(m\ Tj 



(2) 



where the functions $(n) ~ L D (r\n) (n\r') ~ 1. From 
this, one can extract the generic behavior in power 
laws L~ D , which only depends on the space dimen- 
sion in which the system is embedded. The other 
part accounts for phase interferences and for instance 
in the case of almost free electrons (weak disorder 
limit), $(n) = exp(ik.(ri — rj)) so that 5Re$(n)$*(m) ~ 
e -ifc p .(r < -r J ) e +ifcp.(r < -r J ) „ cos(2fe F |r; - r 3 \) given that only 
states (m,n) very close to the Fermi surface will con- 
tribute to the electronic response. Possible limitations of 
this treatment due to strong correlations between elec- 
tronic levels will be discussed later on (see part V.C). 

From our argument, two different effects influencing 
the physical properties may be drawn. On the one hand, 
close to metal-insulator transition peculiar electronic lo- 
calization (like multifractality of eigenstates) may lead to 
correction of the law l/|r, — rj\ D . Indeed, in quasiperi- 
odic systems, quasicrystals or disordered fractal struc- 
tures, the average displacement in time of an electron 
is driven by anomalous diffusion. One can wonders how 
electronic susceptibility is affected by such localization 
effects. 

Besides, if the LDoS is inhomogeneous, i.e if it 
presents strong local fluctuations, an effect of local envi- 
ronment may be deduced from the equation (2) and lead 
to specific features for RKKY coupling. Actually, since 
the LDoS | r) can be expanded in the eigenstates basis of 
H as J2 n $(E ~ e ") I ( r I n ) I 2 ' a substantial increasing of 
the LDoS around the magnetic sites, enhancing the cor- 
responding amplitudes (r \ ri), should also qualitatively 
lead to an increase of the non-local susceptibility. This 
has been discussed in relation with peculiar properties of 
AlPdMn quasicrystalline phases where the RKKY inter- 
action was found to be anomalously large and identical 
for several phases with different TDoS |pd|-[l3|] . 



III. REAL-SPACE APPROACH OF RKKY 
INTERACTION 



with E n > Ef > E m and (n\rj) ~ l/VL D given that 
states are normalized. Given that the average spac- 
ing between energy levels for a given length is roughly 
AE = (piEp)^)' 1 , with p(E) the total density of states, 
one can assume that the energies of the electron and hole 
are respectively E„ — Ef + nAE and E m = Ef — mAE 
(n,m > 0). Consequently, the susceptibility can be ex- 
pressed as 



X (r l ,r J )=2(L- D ) 2 5Re^ 
_ MEf] 
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To perform real space calculations of the non-local 
susceptibility, considering the Green operator G(z) — [z — 
H)^ 1 = J_ °° s( -fS^ dE, one starts from the general form 
oix(n,rj) : 



2Ke / dEdE 

E>E F 



,(n \8{E-H)\r j ){r j \5{E' -H)\n) 
E-E< 



The aim of the method is to determine the coeffi- 
cients (ri | 5(E — H) | Tj) without exact diagonaliza- 
tion, usually limited to simple models and small finite 



size systems. The keypoint of the method is to use a ba- 
sis of orthogonal polynomials {V n (E)} n eN associated to 
a normalized function p(E), referred as a model density 
of states. If the spectral subset of p(E) contains the one 
of the real hamiltonian, it can be shown that 113] 

S(E -H) = P (E)J2'Pn(E)V n (n) 

n 

and the {V n {E)}n€N satisfy to the orthogonality condi- 
tion 



-Hoc. 



P {E)V n (E)V m (E)dE = S„ 



with V-i{H) = and V (H) 
weight is given by 



1 and the corresponding 



Pab(E) 



1 



2vr6 2 
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which is 7^ only for E £ [a— 2b, a+2b\. The coefficients a 
and b are given by the calculated limits a n ~>oo, b n ^oo for 
the real densities of states. From the abovementionned 
relations, the | (pf) will be given by | ipf) = ^(H — a) | 
<)+ | ipr 1 ) and | if' 1 ) = 0,| tf) =| n). In the 
tight-binding scheme H — ^ 7 p9 | r p )(r q |, one shows 
that the a™ coefficients have to be evaluated recursively 



through (a 



= 0,a p] =S pj , 



Vp) 



These relations enable to write {n \ 5(E — H) \ rj) 



p( E )E n a N ^( E X -with a?. 



Vn{H) 



rj ) . From 



these expressions, the susceptibility can be written down 



X(n,rj) = 5Re^Z mn a™a™ 



T 

- L "Tt 



e>Ef p { E)p { E>) V ^f\ E d Ei 

E' <E F 



Accordingly, the calculation of x{ r i: r j) is divided 
into two independent parts. Depending on the choice of 
polynomials, the T mn will display particular analytical 
form. The other part of the susceptibility will involve 
a recursive evaluation of coefficients a?j by means of the 
three term relations defining the orthogonal polynomials, 
as detailed below. 



A. Recursion evaluation of a?j 



One has to choose a convenient basis of orthogonal 
polynomials to evaluate step by step the coefficients afj. 
Given that any basis of such polynomials is uniquely de- 
fine by a three-term recurrence relation, generically read- 
ing EP n (E) = a n V n (E) + b n T n+1 (E) + b n -xP n -x{E) 
with 6_i = 0,n > 0, and a n ,b n the associated recursion 
coefficients Q, the main vectors to be evaluated recur- 
sively will follow from 



Vi 



Pn(H) 



£<# 



Pratically, we will consider the Chebyshev polynomi- 
als of second order that have been already used in others 
contexts. Such polynomials are defined by 

V n (H) = Un - a)Vn{n)+V n -i{H) 






Ipj 



B. Calculation of T mn for the Chebyshev polynomials 



Let us now proceed to the calculation of T mn for 
the Chebyshev polynomials of second order. First we 
rearrange the general form 



P (E)p(E') VmiE) ' PniE ' ) dEdE' 
e>e f PV " y ' e - E' 



by noticing that the factor 1/(2? — E') can be written as 

- e(E - E F )Q{E F - E')) 

with the Heaviside function (0(x) = 0, x < and <d(x) = 
1, x > 0) and the contour T is shown in the complex plane 
on Fig.l for rj — > 0, Rr — » oo. One then rewrites : 



jiEy-Lrtin ■+- -L-nri 
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dz 
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,p( E')V m (E') 
i-E' 



By application of Jordan Lemme, the integral on the 
contour T tends to zero when the radius goes to infinity, 
and it only remains four integrals on the real axis respec- 
tively for [-oo, a -2b], [a-2b,E F ], [E F ,a+2b], [a+26,+oo]. 
Using the relation between first and second order Cheby- 
shev polynomials defined on [—1, +1] and associated with 
p(E) = ^/l^W 



n m i +1 viEEM^m dE = 



*0 = 



u) + ii] — E 



w{Q n+1 (uj) =F i7r\/l - w 2 P„(u>)} 

it is easy to show that for | w |< 1, with a; = coscj) then 

Q„(w) = cos n<p 

sin(n + 1)4> 



Pnfc) 



sin< 



and finally 



f + VI- E 2 P n (a + 2bE) , . . 

lim / ^— ^d£ = 7rexp(=Fi(n+lU 



In conclusion, given that the integration outside 
[— 1, 1] leads to pure imaginary terms, the calculation of 
Imn reduces to 



Ke(T mn +l„ m ) = — — / du)sin((m + n + 2)4>) 

2tto 



Ep-a 
| / 2F" 

2^6 



du; sin((m + n + 2)0) 



1 /sin(m + n + 3)Af sm(rn + n + 1)Af 
2nb \ m + n + 3 m + n + 1 

where Af =Arcos( ' E ^ b ~ a ). The final form of the elec- 
tronic susceptibility for a general tight-binding hamilto- 
nian will be defined by 



1 v~~* n m { sm(m + n + 3)Af sin(m + n + 1)Af 
2irf> 2-*i QiJ aji \ m + n + 3 m + n + 1 



which is the final form of the algorithm. 



IV. DISCUSSION ON THE NUMERICAL 
CONVERGENCE AND TESTS 



This can be done by noticing that n-recursion steps gives 
^ (W the total bandwith) as a resolution in energy, and 
that the states that will mainly contribute to RKKY, 
are enclosed in an interval of width AE < H/t^ around 
Fermi energy. Here, we define tl as the time needed 
for electrons to travel between two magnetic impurities 
distant of L. For a sufficient resolution of our spectrum, 
the number of recursion steps to be considered is n > 
Wtl/1%. 

This result may also be recovered by looking at the 
propagation of recursion states in real space. In fact, 
the propagation of the n-th recursion state turns out 
to be representative of the real wave-function at a time 
t - nh/W, | ipf) -| *(£„ - nh/W)), initially located 
at the same site | r,). Then, one has to check that the 
diffusion length of | </?"), which we define as £(n), is suf- 
ficiently larger than the distance between magnetic sites. 
Therefore as soon as £(n) S> L, the coefficients a^f~ n will 
not contribute significantly. 

The criterion n 3> Wtl/Ti will depend on the nature 
of propagation through the scattering time. For ballistic 
regimes (crystals), tl ~ Vp -L with vp the Fermi energy, 
and n 3> (W/hvp).L, whereas diffusive regimes will be 
associated to L 2 = D.tl, with D the diffusion constant 
and n ^> (W/hD)L 2 . Anomalous diffusion described by 
L ~ At 13 (A,0 < (3 < 0.5, or 0.5 < (3 < I depending on 
the model) will lead ton> (W/hA).]^ 1 ^ to be control 
ipmpirically 

' To test the efficiency of the method, we have checked 
that the susceptibility for tight-binding electrons with 
Fermi energy close to the band edges was equivalent 
to that of free electrons, i.e. well described by the 
law x( r ) ~ (2k F rsm(2k F r) - cos(2fc F r))/(fcFr) 3 (in 2- 
dimension). Random systems have also been simulated, 
by considering site energies distributed at random within 
[— W, W]. For distance larger than the mean free path, 
the decrease of the averaged susceptibility was found to 
be proportional to -^ exp(— r/l pm ), and the numerical 
value of the mean free path was in good agreement with 
the expected one. Calculations were performed for sev- 
eral values of W with average over 150 configurations. 



To achieve the convergence of the calculation, one 
must ensure first of all, that the total length of the system 
is much larger that the distance between magnetic sites 
if one wants to consider an infinite medium. To close the 
system, we take periodic boundary conditions. Second, 
the numerical accuracy must be checked. We found that 
single precision was sufficient and gave the same results 
as double precision. 

It remains to figure out how many recursion steps are 
needed (summation on n and m in the above formula). 



V. APPLICATIONS 

A. 2-D Tight-Binding simulation of a Hume Rothery 
Band Alloy 



In the following, we exemplify the role of local envi- 
ronment on RKKY interaction in a 2D system. Starting 
from a perfect crystal, we consider how perturbation of 
local order in the vicinity of the magnetic site may affect 
RKKY coupling on larger (mesoscopic) scales. The local 



disorder will be designed by a proper modification of lo- 
cal tight-binding parameters of the corresponding hamil- 
tonian. By doing so, we ensure that the LDoS around 
the initial state | ipo) is strengthened in regards to TDoS 
at the same energy. 



Uj> iy) 

+ (% + ij i +<* - ij i )+t J2 1 #>«{* ± y ± 1} 



u> 



Our model (Fig. 2) features two different hopping integral 
between first (7) and second (r) nearest neighbors as well 
as an alternate distribution of site energies 



= t" + ' 



\i+3 



+HL(1_ (_!)<+•' 



in the r = case, this leads to a two-band DOS with a 
real gap whose bandwith A = eb — £a (e.4 < cb)- 

The orthonormal recursion basis \tj)„) and coefficient 
can be evaluated iteratively. Nonetheless, even for simple 
tight-binding models, the complexity rapidly increases 
§J. Taking e A = 0.1, e B = 0.3,7 = 0.1, r = 0.09 has 
TB parameters one gets a TDoS and a RKKY interac- 
tion as shown in the inset of Fig. 3, and one notices in 
particular that the RKKY coupling shows a structure 
that manifests the two underlying periods of the lattice. 

B. Effect of local environment on RKKY 



complex on mesoscopic scales. Assuming that local en- 
vironments of magnetic sites are associated with strong 
local densities of states (when compared to the average 
total density of states), then the RKKY coupling could 
be anomalously strong when compared with other metal- 
lic disordered phases Q. One notes that such strong 
fluctuations of LDoS at a given energy from site to site 
are also thought to be important at the proximity of a 
metal-insulator transition ||. 



C. Anomalous diffusion and RKKY in quasiperiodic 

systems 



General properties of electronic susceptibility in 
quasiperiodic systems are difficult to describe when com- 
pared with periodic ones. Indeed, the simple law ob- 
tained for free itinerant electrons may be questionable 
due to intrinsic incommensurability. We show in this 
section, that the use of recursion method gives here some 
unique informations about RKKY in quasiperiodic struc- 
tures, and further provide a framework for investigating 
mesoscopic interaction in aperiodic systems. Hereafter 
we consider a 2D-Fibonacci quasilattice for which we will 
analyse the susceptibility in relation with the spectral 
properties on larger systems, by varying the length of our 
systems (25.000 sites), number of recursion steps (up to 
500) and the intensity of quasiperiodic potentials. Corre- 
sponding spectral structure may be found elsewhere ]L8| . 
The hamiltonian is defined for a Fibonacci-2D quasilat- 
tice and written in tight-binding basis as 



Let us now elaborate on the nature of the local per- 
turbation around a given magnetic sites. The idea is to 
particularize the local density of states LDoS around one 
given site and to measure the effect on the susceptibility 
on mesoscopic scale. 

Concretely, we modify the fourth first shells around 
the initial site | tJjq) in the following way 7' = 0.65 7, r' = 
1.26 r and the site energies respectively s'a b = £ a,b + 
Xi, i = {1,2,3,4} with Ai = 0.15, A 2 = 0.12, A3 = 
0.09, A4 = 0.06. Initial site energy is ea = (V'ol^lV'o)- 
Thereby, the corresponding LDoS (pi(E) — — 2/-7rQ ! m(ri | 
G(z) I Ti)) is increased if Fermi energy is Ep = 0.176 eV 
(pi(E) > PBuik{E)) as shown on Fig. 4. The number of 
sites considered for the computation of the susceptibility 
in the numerical was about 250.000, and the number of 
recursion steps around ~ 200. 

The calculated electronic susceptibility for the ho- 
mogeneous and inhomogeneous cases are plotted on Fig. 5 
from the results, one clearly sees that such enhancement 
of LDoS leads to an increase of the susceptibility. This ef- 
fect may be at the origin of peculiar magnetic properties 
pd| in quasicrystals where atomic order is known to be 



U = ^£ l \ Vi){Ti 



-7E 



*i)( r j 



e iy , with 



for which site energies are given by e, = Si 
£ix,£iy = ^V qp (potential strength) according to a Fi- 
bonacci sequence. Hopping integrals are set constant for 
simplicity ((ij) denotes first neighbors). On Fig. 6, the 
TDoS for a 2D quasiperiodic Fibonacci quasilattice, as 
well as typical signature of incommensurate long range 
order are reported. The strength of the quasiperiodic po- 
tential is V qp — 0.7(7) an d the susceptibility is given as 
a funtion of the distance between interacting magnetic- 
site (in a- units, with a the lattice spacing). From an 
analysis as a function of the potential strength, we found 
that no Fermi wavelength can be properly defined and 
oscillations exhibit resurgences that are absent from the 
periodic potential. This is a surprising pattern absent 
for periodic potentials for which unique wave-vector (at 
Fermi level) and continuous decreasing of the coupling is 
found. Such patterns may however remember the local 
fluctuations found that we found in random systems for a 
given configuration of disorder. This may thus appear as 



a common feature between the RKKY in quasiperiodic 
and disordered systems on mesoscopic scales. 

Incommensurability effects on the period can be re- 
vealed by analysis of the fluctuations of log \xi r i — r j)\ 
as a function of log(ri — rj). Indeed, if the V qp = 0.0 
case manifests the oscillating behavior cos(2fcp|r|) with 
unique wave-vector, quasiperiodic potential breaks this 
pattern even for small values such as V qp = 0.05(7). O n 
Fig. 7, we compare three results obtain at Ep = — 1.9"f 
and V q p = 0.0,0.05,0.1. Appearance of site-dependent 
incommensurate phase shift is thus illustrated and un- 
veil the action of quasiperiodic potential on electronic 
coherence (log \x\ has been rescaled for more clarity) 

We now consider the long range properties of RKKY 
in these systems. To that end, let us return to the rela- 
tion between electronic localization and diffusion modes 
at the origin of interchange coupling. In the classi- 
cal case (Brownian motion), assuming p(r,t) the prob- 
ability density of finding a random walker (or a classi- 
cal electron) at r after time t, the anomalous diffusion 
regime (found in a disordered fractal) is expressed by 
(r 2 (t)) dta = fr 2 p(r,t)d D r ~ t 2v |]|. Anomalous regime 
can also be found in quantum diffusion, especially in 
quasiperiodic systems |9[jl0|] and at the metal-insulator 
transition (quantum Hall systems) fllq] . This has been 
connected with multifractal properties of eigenstates. In- 
deed if we consider the spreading of wave-packet con- 
structed from multifractal eigenstates, one finds numer- 
ically that (f 2 ) = <*(i)|r 2 |*(i)} = f r 2 \1>(r,t)\ 2 d D r ~ t 2v 
with v an exponent characterizing the strength of the 
potential, the intrinsic correlations, etc. According to 
an argument by E. Akkermans |l7| ] the power-law de- 
creasing of the averaged second moment of electronic 
susceptibility should not be affected by anomalous dif- 
fusion and remains of type (x(\( r i — r j\)) ~ \ r i — r j\~ 
for D-dimensional disordered systems. 

For the quasiperiodic strengthes considered in our 
calculation (from V qp = 0.05(7) U P to V qp = 1.1(7)), the 
long-range power law is not affected qualitatively and 
one recovers the l/|r|~-° with D = 2. This is exempli- 
fied on Fig. 8 for two different intensities V qp = 0.1 and 
Vqp = 1.1 compared with the long ranged oscillations for 
no quasiperiodic potential which mimic the periodic po- 
tential. Note that the slight departure from power-law 
in the case V qp = 1.1 (and for large distance) is due to 
finite-size effects and can be smoothed out by increasing 
recursion steps. 

One remarks that results on the octagonal quasiperi- 
odic tiling have been obtained previously p9| , and also 
seem to indicate a significant site dependence of the elec- 
tronic susceptibility apparently without alteration of the 
general power-law as a function of distance between mag- 
netic sites. 



These results are in agreement with the argument 
we gave in section II for the general dependence of 
X (n, rj ) = (2p(E F ))/L D ^eJ2 n , m HnW(m)/(n + m) which 
gives a correct interpretation of the universal shape of 
L~ D for a metallic system. But our assumption made 
on the level distribution does not include any subtle cor- 
relations in the spectrum, as those found at the metal- 
insulator transition and associated to a specific level spac- 
ing distribution EQ] . If quasiperiodic potentials induces 
multifractal states, their general level spacing distribu- 
tion has been however found to be described by the 
Gaussian Orthogonal Ensemble (GOE) of the Random 
Matrix theory pl}| . This involves a metallic state of 
the conduction system and our calculations are conse- 
quently in agreement with other recent results obtained 
on quasiperiodic systems. Manifestations of multiffac- 
tality of eigenstates may be rather revealed by strong 
fluctuations of local densities of states || . 

To finish with, one notes that unexpected feature 
has been also revealed from careful analysis of site de- 
pendent susceptibility (to be published elsewhere). For 
a potential V qp = 0.1 and Fermi energies Ep = —0.65(7) 
(resp. Ep = —0.575(7)) situations where the electronic 
coupling is purely ferromagnetic (resp. antiferromag- 
netic) were found p2[ . Such kind of pattern unveil un- 
precedented signature of complicated localization effects 
unique to quasiperiodic structures. Related phenomena 
such as Kondo effect in quasiperiodic systems has been 
discuss in regards to the same localization effects p2| . 



VI. CONCLUSION 



A method to investigate RKKY interaction in pe- 
riodic or aperiodic systems has been presented. Using 
real-space schemes, this approach enables in particular to 
analyse the effect of local inhomogeneities, quasiperiodic 
potential and by extension in all situations where usual 
diagonalization methods may be limited. Some conclu- 
sive check of the method have been given in comparison 
to expected behavior and specific models have been stud- 
ied in order to exemplify the interest of our method. In 
particular we have shown that even strong regimes of 
localization induced by a quasiperiodic potential do not 
lead to a qualitative departure from general power-law 
dependence, and LDoS fluctuations could affect the in- 
tensity of the interaction. 
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VIII. FIGURES CAPTIONS 



Fig. 1. Integral contour for the calculation of the coeffi- 
cient Iran and semi-elliptic density of states p a b(w) used 
in the case of Chebyshev polynomials. 
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Fig. 2. Schematic representation of the lattice with dif- 
ferent tight-binding parameters. 

Fig. 3. Total density of states for the hamiltonian in the 
perfect Hume-Rothery alloy and corresponding RKKY 
coupling (inset). 

Fig. 4. Local density of states around the modified local 
environment of one magnetic site. 

Fig. 5. Comparison of the electronic susceptibility (in 
a- units, a lattice parameter) for the homogeneous (bold 
line) and the inhomogeneous cases. 

Fig. 6. RKKY coupling for a quasiperiodic system 
(V qp = 0.7) and different Fermi energies. The inset shows 
the corresponding total density of states. 

Fig. 7. Electronic susceptibility as a function of 
log |x| for different small quasiperiodic potentials (V qp = 
0,0.050.01) 

Fig. 8. (a) Electronic susceptibility as a function 
of log |x| for small and strong quasiperiodic potentials 
(Vqp = 0.11.1) compared to periodic potential. The in- 
set (b) shows the local susceptibility for V qp = 0.1, the 
same Fermi energy and for several different environment 
between magnetic sites in the quasiperiodic potential. 
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